The Phosphoproteome of the Rd1 Mouse Retina, a Model of Inherited Photoreceptor Degeneration, Changes after Protein Kinase G Inhibition

Retinitis pigmentosa (RP) is a frequent cause of blindness among the working population in industrial countries due to the inheritable death of photoreceptors. Though gene therapy was recently approved for mutations in the RPE65 gene, there is in general no effective treatment presently. Previously, abnormally high levels of cGMP and overactivation of its dependent protein kinase (PKG) have been suggested as causative for the fatal effects on photoreceptors, making it meaningful to explore the cGMP-PKG downstream signaling for more pathological insights and novel therapeutic target development purposes. Here, we manipulated the cGMP-PKG system in degenerating retinas from the rd1 mouse model pharmacologically via adding a PKG inhibitory cGMP-analogue to organotypic retinal explant cultures. A combination of phosphorylated peptide enrichment and mass spectrometry was then applied to study the cGMP-PKG-dependent phosphoproteome. We identified a host of novel potential cGMP-PKG downstream substrates and related kinases using this approach and selected the RAF1 protein, which may act as both a substrate and a kinase, for further validation. This showed that the RAS/RAF1/MAPK/ERK pathway may be involved in retinal degeneration in a yet unclarified mechanism, thus deserving further investigation in the future.


Introduction
The hereditary disease Retinitis pigmentosa (RP) is a frequent cause of blindness among the working population in industrial countries [1]. There is no effective treatment available at present [2], with the notable exception of the approved gene therapy for mutations in the RPE65 gene [3]. RP results in vision loss via fatalities of the light-sensing photoreceptors, which are categorized as rods and cones. The degeneration pattern starts with the death of rods that cause night blindness, which is followed by the degeneration of cones [4] and hence a loss of general vision. The rd1 mouse model, with the gene mutation of the beta subunit of the enzyme phosphodiesterase-6, is characterized by an abnormally high level of photoreceptor cyclic guanosine monophosphate (cGMP) and a rapid loss of these cells [5]. The degeneration of photoreceptors in this phenotype has been demonstrated with an increasing decline in outer retinal thickness and overall amplitude reduction in the waveform via spectral domain optical coherence tomography [6] and electrophysiological changes [7], respectively. Although the details on how the high level of cGMP affects the photoreceptors are still not known, previous studies have suggested that at least part of the death-promoting impact on photoreceptors could be exerted via increased activity of cGMP-dependent protein kinase (PKG) probably causing over-phosphorylation within photoreceptors [8,9]. It would thus be most helpful if this could be studied more closely.

Quantitative Phosphoproteomics of Retinal Explants with PKG Inhibition
In total, 17,915 phosphorylated sites and 3405 retinal proteins were identified ( Figure 1A); among these, 10,352 phosphosites from 2214 proteins were selected for further analysis, according to the defined criteria (detailed in Section 4). Comparison via Student's t-test showed that 299 phosphosites of 108 proteins had decreased phosphorylation, while 184 phosphosites of 83 proteins had increased phosphorylation in retinas with PKG inhibition for 2 h compared to their untreated counterparts ( Figure 1B, Table 1; see Supplementary  Table S1 for a full list, with phosphorylation sites, fold changes, etc.). For the longer treatment, 4 days, 98 phosphosites of 46 proteins were observed with decreased phosphorylation, while 265 phosphorylation sites of 130 proteins with increased phosphorylation. ( Figure 1C, Table 2, and Supplementary Table S2 for full list).
There were five overlapping proteins with decreased phosphorylation after both 2 h (short term) or 4 days (long term) PKG inhibition. ( Figure 1D and Table 1; overlapping proteins highlighted in yellow in Supplementary Tables S1 and S2). For proteins with increased phosphorylation, there were 16 overlapping proteins after short and long term PKG inhibition ( Figure 1E and Table 2; overlapping proteins highlighted in yellow in Supplementary Tables S1 and S2). When we focused on individual phosphorylation sites in the peptides of the overlapping proteins, there were two proteins with decreased phosphorylation that were represented by the very same, specific phosphorylation sites at both time points, and three proteins from those with increased phosphorylations (Tables 1 and 2; plus labeled with green in Supplementary Tables S1 and S2). Table 1. Lists of proteins showing decreased phosphorylation after PKG inhibition for either 2 h or 4 h, as well as reduced phosphorylation at both time points. The left and center columns only show the 20 proteins with the greatest fold change for 2 h and 4 days, respectively. The total number of proteins with decreased phosphorylation was 108 and 46 after 2 h and 4 days of PKG inhibition, respectively. Full lists of all proteins/peptides with decreased phosphorylation are found in Supplementary Tables S1 and S2. phosphorylation that were represented by the very same, specific phosphorylation sites at both time points, and three proteins from those with increased phosphorylations (Tables 1 and 2; plus labeled with green in Supplementary Tables S1 and S2).

Potential cGMP-PKG Dependent Biological Pathways
A protein annotation in the Kyoto Encyclopedia of Genes and Genomes (KEGG), with the targets mentioned above, was performed to address the biological pathways in which the cGMP-PKG-dependent substrates were involved. Proteins with decreased phosphorylated sites after PKG inhibition were annotated in RNA transport, Spliceosome, and Rap1 signaling pathway ( Figure 2A). A part of these down-regulated sites after PKG inhibition, such as the sites in protein kinase cGMP-dependent 1 (PKG1), general transcription factor IIi (GTF2I), and protein phosphatase 1 catalytic subunit alpha (PPP1CA), were involved in the cGMP-PKG pathway. This supported that PKG signaling in explants was appropriately inhibited during culturing. By contrast, the proteins of the increased phosphorylation were mostly related to pathways connected with proteoglycans, ErbB signaling, and tight junctions ( Figure 2B). In this manner, we identified potential signal pathways that may function under cGMP-PKG regulation in photoreceptors undergoing degeneration.

Upstream Kinase Analysis
To predict the upstream kinases for the altered phosphorylations identified above, the regulated phosphoproteins were analyzed via annotation to the online tool Phosphomatics V2 (https://phosphomatics.com/). The outcomes are given as kinome trees ( Figure  3A,B) and as lists (Supplementary Tables S3 and S4). For the 2 h treatment, the algorithms suggested a total of 66 potential upstream kinases, 32 with decreased activity and 34 with increased activity ( Figure 3A, Supplementary Table S3). Among these, 7 kinases in the category of proposed decreased activity, and 11 of the increased ones phosphorylate their substrates in a site-specific manner (Supplementary Table S3).
For the 4 days treatment, the number of suggested upstream kinases reached 68, of which 8 had decreased activity and 60 had increased activity ( Figure 3B, Supplementary  Table S4). In this case, the site-specific ones were represented by three and five kinases of those with proposed decreased and increased activity, respectively (Supplementary Table  S4). The kinases that had a consistent reaction to the treatment between the two time points were summarized (Table 3). Two kinases gave indications of having reduced activity after both time points, namely cyclin-dependent kinase 1 (CDK1) and cAMP-dependent protein kinase catalytic subunit alpha (PRKACA), whereas 10 kinases came up as having increased activity at both time points. Of the latter, four kinases were related to the mitogen-activated protein kinase (MAPK) family (Table 3). There were also some kinases that displayed opposite reactions to the treatment after 2 h and 4 days (Table 3), such as CDK2, CDK7, MAPK3, and more. These are not considered further in this report. To sum up, according to the upstream kinase analysis with the selective differentiated

Upstream Kinase Analysis
To predict the upstream kinases for the altered phosphorylations identified above, the regulated phosphoproteins were analyzed via annotation to the online tool Phosphomatics V2 (https://phosphomatics.com/). The outcomes are given as kinome trees ( Figure 3A,B) and as lists (Supplementary Tables S3 and S4). For the 2 h treatment, the algorithms suggested a total of 66 potential upstream kinases, 32 with decreased activity and 34 with increased activity ( Figure 3A, Supplementary Table S3). Among these, 7 kinases in the category of proposed decreased activity, and 11 of the increased ones phosphorylate their substrates in a site-specific manner (Supplementary Table S3).
For the 4 days treatment, the number of suggested upstream kinases reached 68, of which 8 had decreased activity and 60 had increased activity ( Figure 3B, Supplementary  Table S4). In this case, the site-specific ones were represented by three and five kinases of those with proposed decreased and increased activity, respectively (Supplementary Table S4). The kinases that had a consistent reaction to the treatment between the two time points were summarized (Table 3). Two kinases gave indications of having reduced activity after both time points, namely cyclin-dependent kinase 1 (CDK1) and cAMPdependent protein kinase catalytic subunit alpha (PRKACA), whereas 10 kinases came up as having increased activity at both time points. Of the latter, four kinases were related to the mitogen-activated protein kinase (MAPK) family (Table 3). There were also some kinases that displayed opposite reactions to the treatment after 2 h and 4 days (Table 3), such as CDK2, CDK7, MAPK3, and more. These are not considered further in this report. To sum up, according to the upstream kinase analysis with the selective differentiated phosphoproteins mentioned above, a series of kinases regulated by cGMP-PKG were revealed to potentially function during retinal degeneration.

Selection of RAF-1 for Further Validation
To gain more insights in the cGMP-PKG network, one target was selected from the upstream kinase analysis for further validation. Since the abnormally high level of cGMP in the photoreceptors is regarded as a disease driver of retinal degeneration [19], we assumed that photoreceptor-specific proteins have an essential role in retinal dysfunction.

Selection of RAF-1 for Further Validation
To gain more insights in the cGMP-PKG network, one target was selected from the upstream kinase analysis for further validation. Since the abnormally high level of cGMP in the photoreceptors is regarded as a disease driver of retinal degeneration [19], we assumed that photoreceptor-specific proteins have an essential role in retinal dysfunction. For this purpose, the RAF proto-oncogene serine/threonine-protein kinase (RAF1) was selected for additional analysis. As mentioned above, several of the upstream kinases with increased activities belonged to the MAPK family. The rationale for selecting RAF is thus that MAPK1 is expressed in photoreceptors and is linked to photoreceptor death, at any rate in Drosophila [20], and also that RAS and its downstream pathway component RAF1 are required to satisfy photoreceptor differentiation [21]. In addition, RAF1 was proposed as having higher kinase activity after 2 h of PKG inhibition (Supplementary  Table S3). Therefore, as part of the well-established Ras/Raf-1/MAPK1 signaling, the role of RAF1 during retinal degeneration and its connection with the cGMP-PKG system was of interest to analyze further. Other kinases identified in this study appear to be good subjects for deeper investigation in the future.
To see whether RAF1 is specifically expressed within photoreceptors, retinas from either rd1 or wt animals were immunostained for RAF1. In wt, RAF1 was expressed predominantly in the photoreceptor segments. This was also true for rd1, which in addition, displayed elevated RAF1 staining in many parts of the rd1 outer nuclear layer (ONL), where the photoreceptor cell bodies and processes reside ( Figure 4A-C). This supports that RAF1 can be seen as a relatively photoreceptor-specific kinase and that it could be involved in photoreceptor degeneration.
Next, the potential link between RAF1 expression and the cGMP-PKG system was investigated at the histological level. This was achieved via pharmacological PKG inhibition in retinal explants during culturing (the same culturing paradigms as above), followed by an analysis of possible RAF1 alterations in photoreceptors after treatment. Figure 5 shows that the RAF1 expression was higher in the ONL after PKG inhibition, regardless of the different lengths of treatment.   To see if the change in RAF1 expression was in any way also reflected in the phosphorylation state of RAF1, the PLA technique for the combination of RAF1 protein and phosphorylated serine [22] to study the RAF1 phosphorylation profile was applied. The PLA uses antibodies against two different antigens, as well as certain secondary antibodies containing complimentary "primers" of a reaction product. If the two antigens are in close proximity (<40 nm), the method then allows a reaction product to be formed. Thus, if the RAF antibody is in close contact with the phosphoserine antibody, it is likely that both antibodies sit on the same protein, namely RAF1, and thus identify phosphorylated RAF1. We focused on the RAF1 phosphorylation in ONL since RAF1 was higher expressed in the rd1 ONL ( Figure 4) and found an increased level of PLA signal within the ONL of rd1 ( Figure 6). This indicated a greater extent of RAF1 phosphorylation in retinas undergoing degeneration.  The next question was if the RAF1 phosphorylation could be affected by the manipulation of the cGMP-PKG system during culturing. We observed that RAF1 phosphorylation was higher in the ONL of retinal explants being treated with PKG inhibitor for 2 h ( Figure 7A,B,E), while no difference to the untreated counterpart was noticed after 4 days of treatment ( Figure 7C-E). This was consistent with our discovery in the upstream kinase analysis that a higher RAF1 activity was identified in explants with 2 h of PKG inhibition but not in the 4 days treatment ( Figure 3A,B and Supplementary Tables S3 and S4). The combined data thus clearly indicate the participation of RAF1 activities in retinal degeneration during the 2 h of PKG inhibition, albeit in an unclarified mechanism. It is reasonable to deduce that different lengths of PKG inhibitory treatment affect the RAF1 activities, such that RAF may be only transiently affected by the change in PKG activity, but that this still results in effects downstream to RAF1 itself. Note that the overwhelming amount of the phosphorylated RAF1 resided in the ONL and the photoreceptor segments, again suggesting a potential role of RAF1 during retinal degeneration, though further investigation is required to specify whether RAF1 phosphorylation promotes neurodegeneration or neuroprotection in photoreceptors.
= 4) and their untreated controls (n = 3). (A,B) represent RAF1 (green) immunostaining in untreated retinal explants and their counterparts with 2 h of PKG inhibition. (C,D) represent immunostaining of RAF1 in untreated retinal explants and their counterparts with 4 days of PKG inhibition. DAPI (blue) was used as nuclear counterstain. (E) The box chart compares RAF fluorescence intensities within the ONL between untreated explants and peers with PKG inhibitors in different treatment lengths. * p < 0.05. Figure 6. Comparison of RAF1 phosphorylation between rd1 and wt retinas. The PLA punctuations indicate that RAF1 is in proximity with phosphoserine, which reveals the RAF1 phosphorylation profile. (A,B) represent PLA punctuations (RAF1 and phosphorylated serine) in ONL from rd1 and wt, respectively. DAPI (blue) was used as nuclear counterstain (C) The box chart shows the comparison of PLA counts in ONL of retinas from rd1 (n = 5) and wt (n = 5). **** p < 0.0001.

Discussion
Despite gene therapy targeting the mutated RPE65 gene having been developed to treat RP [3], other effective treatments are required because of the genetic heterogeneity of the disease [23]. The cGMP-PKG system [5], which is regarded as one of the RP disease drivers, is the potential molecular target for a novel therapeutic strategy development for this disease. In this context, the cGMP-PKG system is supposed to function via the phosphorylation of retinal substrates, including in photoreceptors [8,9]. It is, therefore, likely that cGMP-PKG-dependent phosphoproteomics will provide more insights into the details of retinal degeneration. In this case, MS-based phosphoproteomics techniques, with extensive application in biological research [24], would be a proper method to investigate the cGMP-PKG downstream signaling. Technically, one option is to apply immunoprecipitation for phosphoprotein enrichment prior to MS detection. However, due to the limitation of commercially available phospho-selective antibodies, this approach has not provided promising results in phosphoprotein enrichment [25] and is, therefore, unsuitable for high throughput phosphoproteomic study. Another widely applied strategy is the use of immobilized metal affinity chromatography, as in the current study. Though the non-specificity has been noticed for this method, the continuous optimization of the protocol has made it frequently applied in phosphoproteomic research [26], including here.
Our study has three clear outcomes: (1) it identified potential retinal PKG substrates that could be involved in retinal degeneration, (2) it gave information on the kinase and pathway networks in operation during the degenerative events, and (3) it showed a correlation between the cGMP-PKG system and the RAS/RAF/MAPK/ERK pathway in this situation.
The upstream kinase and KEGG analyses both picked up PKG or the cGMP-PKG signaling pathway as being affected by the treatment, which corroborated the usefulness of the approach. In addition, the two treatment regimes, i.e., 2 h or 4 days of PKG inhibition, gave data on proteins whose phosphorylation was decreased as well as such that displayed increased phosphorylation in either or both of these regimes.
The assumption that increased phosphorylation by PKG is part of the degeneration mechanism makes the proteins with reduced phosphorylation at both time points of PKG inhibition particularly interesting with respect to a pathologic involvement. Among the five proteins of this category, we noticed DNA Topoisomerase II Beta (TOP2B) and capicua transcriptional repressor (CIC). Interestingly our previous transcriptome study revealed the association between cGMP-PKG and these two targets; hence, there is reason to keep TOP2B and CIC as candidates for links between PKG activity and events during retinal degeneration. Proteins that instead showed persistent increased phosphorylation after PKG inhibition are, by contrast, more likely to be dependent on kinases under direct or indirect negative regulation by PKG and were represented by as many as 16 proteins. Some deserve to be mentioned here, such as ankyrin 2 (ANK2), with a possible role in maintaining ion balance, via anchoring of ion transporters to the photoreceptor cell membrane [27]. Proteins of the microtubule-associated protein (MAP) family, including MAP1B, MAP2, and MAP4, were also seen. MAPs could be critical for photoreceptors, and at least one of these may associate with retinal degeneration [28,29].
The upstream kinase analysis indicated integration of the cGMP-PKG pathway with other kinases. With the same rationale mentioned above, we noticed CDK1 that has a potentially lower activity after 2 h and 4 days of PKG inhibition. Given that we recently demonstrated its expression in degenerating rd1 photoreceptors, as well as that this was reduced by PKG inhibition [30], a possible interpretation would be that CDK1 expression is somehow affected by a cGMP-PKG regulated phosphorylation. By contrast, several kinases with potential higher activities after PKG inhibition were from the MAPK family, which appear connected with the degeneration, including in a PKG-related way.
The RAS/RAF1/MAPK/ERK pathway, with its diverse biological functions, is one of the most extensively studied signal transduction routes [31]. The role of RAF1 has been studied in diabetic retinopathy [32], and here we provide insights into how RAF1 may connect with RP, which has seldom been noticed previously. In wt retinas, the predominant expression of RAF1 was in the photoreceptor segments, pointing to it being a photoreceptor-specific kinase. The increased expression and phosphorylation of RAF1 in the rd1 compared to wt ONL revealed that it may participate in the degeneration of the photoreceptors, not the least since the above suggested increased activation of several MAPKs.
However, neither its increased expression nor the increased phosphorylation specifies whether RAF1 works for or against the degeneration. Our results of PKG inhibition in rd1 explants provided more clues, though, since RAF1 expression increased after both lengths of treatment and more phosphorylated RAF1 was seen after 2 h of PKG inhibition. Given that retinal degeneration benefits from such PKG inhibition [8,9], RAF1 may exert neuroprotective effects during photoreceptor death. Had it been the opposite, that the expression and phosphorylation state of RAF1 drives the degeneration, both these parameters would have been expected to decrease after PKG inhibition. Yet the situation is not straightforward since phosphorylation of RAF1 may lead to either its inhibition or activation, depending on the actual serine site [33], which we were not able to discriminate between here. While highly speculative, the suggested activation of RAF1 in the upstream kinase analysis could mean that the degeneration initiates a protective program involving RAF1 activity but that PKG counteracts the same by blocking the activity, probably via other kinase or even phosphatase activities.
In conclusion, our study shows an intricate pattern of changes with respect to protein phosphorylations after inhibition of PKG. It appears likely that several of these results, including those connected to the RAS/RAF1/MAPK/ERK pathway, relate to the yet unclarified mechanisms of inherited retinal degeneration and, as such, will aid in designing future investigations in this area.

Animals
The C3H/rd1/rd1 (rd1) and control C3H wild-type (wt) were kept and bred in-house under standard white cyclic lighting, with free access to food and water [34]. The animals were used irrespective of sex since it is practically difficult to discriminate at the early age when they were sacrificed in our protocol. The day of birth of the animal was considered as postnatal 0 (P0), with the day following this considered as P1, etc.

Organotypic Retinal Explant Culture
Retinas from P5 mice were used to generate explants according to our standard protocol [17]. Inserts with the explants were put into six-well culture plates with 1.5 mL serum-free medium in each well. Plates were incubated at 37 • C with a 5% CO 2 atmosphere, and the medium was replaced every two days. No additions were made to the cultures for the first 2 days. One group of the cultures, which at this point were at an age corresponding to P7, were then exposed to 50 µM Rp-8-Br-PET-cGMPS (PKG inhibitor; Biolog, Bremen, Germany, Cat. No.: P 007) for the following 4 days, with the endpoint thus equivalent to P11 (n = 4), with their corresponding controls (n = 4) receiving the same amount of distilled water. Another group (n = 4) was exposed to the Rp-8-Br-PET-cGMPS with the same concentration but only for 2 h before the end of the protocol, with controls (n = 4) similarly exposed to distilled water. All retinal explants were either collected for protein extraction, followed by phosphoproteomic analyses. In parallel, another four groups of rd1 explants, reaching a total of 13 explants, were generated. One group of retinas (n = 4) was exposed to the 50 µM Rp-8-Br-PET-cGMPS for 4 days, while their corresponding untreated controls (n = 3) received the same amount of distilled water. The remaining two groups included one (n = 3) that was treated with 2 h of PKG inhibition with the same concentration and one with its untreated peers (n = 3). These explants were fixed, sectioned, and used for microscopy-based studies.
Animals were randomly assigned to the experimental groups, and the experimenters were aware of the conditions of the animals during retinal explant culture. No statistical method was applied to pre-determine the sample size of the experimental groups. The sample size was chosen based on the relevant literature in the field as well as on ethical considerations. The study was not pre-registered. The study was exploratory without exclusion criteria being pre-determined.

Sample Preparation for MS
A total of 16 retinal explants were used for MS measurement, and 4 explants were used per group. The experimenter was unaware of the explant's group during experimentation. Each explant was homogenized separately in buffer (50 mM Tris-HCl, 50 mM NaCl, 1 mM EDTA, 5 mM NaH 2 PO 4 , 1 mM DL-Dithiothreitol (DTT)), supplemented with phosphatase inhibitors (

MS Acquisition and Analysis
The peptide analyses were performed on a Q Exactive HFX mass spectrometer (Thermo Scientific) connected to an EASY-nLC 1200 ultra-high-performance liquid chromatography system (Thermo Scientific). Peptides, 1 µg, were separated on an EASY-Spray column (Thermo Scientific; ID 75 µm × 50 cm, column temperature 45 • C) operated at a constant pressure of 800 bar. A two-step gradient of buffer B (80% acetonitrile, 0.1% formic acid) in buffer A (aqueous 0.1% formic acid) was applied at a flow rate of 300 nL min −1 . In the first step, a gradient of 10 to 30% of buffer B was run for 90 min, followed by a 30 to 45% gradient of buffer B in 20 min. One full MS scan (resolution 60,000 @ 200 m/z; mass range 350-1400 m/z) was followed by MS/MS scans (resolution 15,000 @ 200 m/z). The precursor ions were isolated with 1.3 m/z isolation width and fragmented using higherenergy collisional-induced dissociation at a normalized collision energy of 28. Charge state screening was enabled, and singly charged ions as well as precursors with a charge state above six were rejected. The dynamic exclusion window was set to 10 s. The automatic gain control was set to 3 × 10 6 for MS and 1 × 10 5 for MS/MS, with ion accumulation times of 45 ms and 60 ms, respectively. The intensity threshold for precursor ion selection was set to 1.7 × 10 4 .
The raw DDA data were analyzed with Proteome Discoverer™ Software (Version 2.3, Thermo Fisher Scientific). Peptides were identified using both SEQUEST HT and Mascot against the UniProtKB mouse database (UP000000589 plus isoforms). The search was performed with the following parameters applied: static modification: cysteine carbamidomethylation and dynamic modifications: N-terminal acetylation. Phosphorylation (S, T, Y, for serine, threonine, and tyrosine, respectively) was set as a variable for the phosphopeptide analysis. Precursor tolerance was set to 10 ppm and fragment tolerance to 0.02 ppm. Up to 2 missed cleavages were allowed, and Percolator was used for peptide validation at a q-value of a maximum of 0.05. Peptides with different amino acid sequences or modifications were considered unique. Extracted peptides, with a quality q-value lower than 0.05, and modification of S, T, and Y phosphorylation detected were used to identify and quantify them by label-free relative quantification. The extracted chromatographic intensities were used to compare peptide abundance across samples.
The MS results were processed via Perseus software (version 1.6.0.7) [35]. The protein intensities were log2 transformed. Selective criteria were defined as follows: among 4 replicates of 4 conditions (16 samples total), only a peptide with the missing value of less than 30% in total (less than 5 samples with missing value) would be selected, and the missing values were replaced from a normal distribution was performed through data imputation by using the following settings: width 0.3 and downshift 0. Further bioinformatics analysis of these processed data was performed via the webbased tool Phosphomatics (https://phosphomatics.com/; [36]). In Phosphomatics, two-Sample Student's t-test (two-tailed) was performed to compare phosphorylated site levels between the rd1 explants with 2 h of PKG inhibition vs. untreated control and 4 days of PKG inhibition vs. untreated controls. A p-value of 0.05 was defined as the cutoff. Two lists of the phosphorylated sites identified previously were used to perform a further upstream kinase analysis in Phosphomatics. The biological pathways that might be affected between rd1 untreated and rd1 with PKG inhibition samples were determined in Enrichr (https://maayanlab.cloud/Enrichr/), where differentially phosphorylated sites mentioned above were used as an input list to perform pathway-based analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG).

Cryosection, Immunohistochemistry, and Proximity Ligation Assay (PLA)
Retinal tissues from rd1 and wt in vivo at P9, as well as P11 cultured explants from rd1, were treated with 4% formaldehyde for 2 h, washed 3 × 15 min in phosphate-buffered saline (PBS), cryoprotected in PBS + 10% sucrose for overnight at 4 • C and subsequently with PBS + 25% sucrose for 2 h. After embedding in a medium with 30% bovine serum albumin (BSA; Cat. No.: A5253-250G; Sigma-Aldrich, St. Louis, MO, USA) and 3% gelatine (Cat. No.: 1040781000, Merck Millipore, Burlington, MA, USA) mixed in H 2 O, 12 µm thick retinal cross-sections were cut and collected from an HM560 cryotome (Microm, Walldorf, Germany). The sections were stored at −20 • C for later usage. Cryosections were then used for immunostaining and PLA. Totally 10 retinal samples were used, with 5 each in either rd1 or wt strains, while 13 rd1 explants were used, with 4 in the group of retinas with 4 days of PKG inhibition and 3 each in the other 3 groups.
For immunostaining, briefly, the cryosections were dried at room temperature for 15 min and rehydrated in PBS. Then they were blocked with 1% BSA + 0.25% Triton X100 + 5% goat serum in PBS at room temperature for 45 min. The primary antibody anti-RAF1 (Cat#: MA5-17162, ThermoFisher) was diluted with 1% BSA and 0.25% Triton X100 in PBS (PTX) and incubated at 4 • C overnight; a no primary antibody control ran in parallel. Sections were washed 3 × 5 min each in PTX and incubated with a donkey anti-mouse IgG (H+L) highly cross-adsorbed secondary antibody, Alexa Fluor™ Plus 488 (#A32766, ThermoFisher) at 1:800 dilution in PTX. After 3 × 5-min PBS washes, the sections were mounted with Vectashield DAPI (Vector, Burlingame, CA, USA).
The PLA, which detects if two antigens are in close proximity (<40 nm) with each other, was performed on cryosections, using Detection Reagents Red kit (DUO92008-100RXN, Merck, Readington Township, NJ, USA), the PLA Probe anti-rabbit PLUS (DUO92002, Merck, Readington Township, NJ, USA), and anti-mouse MINUS (DUO92004, Merck, Readington Township, NJ, USA). The procedure followed the manufacturer's instructions. In short, the retinal sections were blocked with blocking solution (provided in PLA Probe kit) for 45 min at 20 • C, and primary antibodies antibody anti-RAF1 (Cat#: MA5-17162, ThermoFisher) and anti-Phosphoserine (ab9332, Abcam, Cambridge, UK) were incubated overnight at 4 • C. PLA Probe anti-rabbit PLUS, and anti-mouse MINUS were incubated for 1 hr at 37 • C. The ligation and amplification steps were performed using the Detection Reagents Red, followed by the addition mounting of Vectashield DAPI (Vector, Burlingame, CA, USA).

Microscopy and Image Processing
A Zeiss Imager Z1 Apotome Microscope (Zeiss, Oberkichen, Germany) with a Zeiss Axiocam digital camera was used for microscopy observations. Image generation and contrast enhancement were performed identically for all images via the ZEN2 software (blue edition). The immunostaining was analyzed for staining differences via three sections, each from three to five animals for each condition, after which the fluorescent intensities of positive cells were randomly distributed within the area of interest (outer nuclear layer, ONL, i.e., the photoreceptor layer) were assessed. Fluorescence intensity was captured and analyzed by the ImageJ software (version 1.53a, NIH, Rockville, MD, USA). The freehand selection function was used to target the ONL, after which the fluorescence intensity was calculated with the measure function. The values of all sections from the same animal were averaged. For PLA, the images were also generated with the aid of a Zeiss Imager Z1 Apotome Microscope and the ZEN2 software (blue edition). The punctuations of ONL, which represent reaction products from when the two antigens are in close proximity, were counted manually in three sections, each from three to five animals for each condition, with values from the same animal averaged. The experimenter was unaware of the groups of the samples during image processing.

Statistical Analysis
The Student's t-test was used to compare the values of immunostaining intensities and PLA punctuations in different conditions in R, and the p-value cutoff was defined as 0.05. Visualization of KEGG analysis was also performed in R. Institutional Review Board Statement: The animal study protocol was approved by the Institutional Review Board (or Ethics Committee) and performed following the issued local animal ethics committee (permit #M92-15 and 02124/2020).

Informed Consent Statement: Not applicable.
Data Availability Statement: The data presented in this study are available on request from the corresponding author.
Conflicts of Interest: P.E. has financial interest in the company Mireca Medicines, which intends to forward clinical testing of the compound used here for PKG inhibition. J.Z. and C.W. have no conflict of interest.